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I discuss the motivation for resummation of the effects of initial-state soft gluon 
radiation, to all orders in the strong coupling strength, for processes in which the 
near-threshold region in the partonic subenergy is important. I summarize the 
method of "perturbative resummation" and its application to the calculation of 
the total cross section for top quark production at hadron colliders. Comments are 
included on the differences between the treatment of subleading logarithmic terms 
in this method and in other approaches. 



1 Introduction and Motivation 

In inclusive hadron interactions at collider energies, tt pair production proceeds 
through partonic hard-scattering processes involving initial-state light quarks q 
and gluons g. In lowest-order perturbative quantum chromodynamics (QCD), 
at 0(a 2 s ), the two partonic subprocesses are q + q — > t + t and g + g — > t + t . 
Calculations of the cross section through next-to- leading order, 0{o? s ), involve 
gluonic radiative corrections to these lawest-order subprocesses as well as con- 
tributions from the q + g initial state p. In this paper, I describe calculations 
that go beyond fixed-orcLet-jjierturbation theory through resummation of the 
effects of gluon radiation aaa to all orders in the strong coupling strength a s . 
The physical cross section is obtained through the factorization theorem 
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a ij (S,m) = —J dri$ij(ri,fji)&ij(Ti,rn,iJ,). (1) 

The square of the total hadronic center-of-mass energy is S, the square of the 
partonic center-of-mass energy is s, m denotes the top mass, /i is the usual 
factorization and renormalization scale, and <&jj (77, fi) is the parton flux. The 
variable 77 = — 1 measures the distance from the partonic threshold. The 
indices ij € {qq, gg} denote the initial parton channel. The partonic xross 
section &ij(rj, m, /_t) is obtained either from fixed-order QCD calculations □, or, 
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as described here, from calculations that include of resummation n un to all 
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orders in a s . I use the notation a = a(fi = m) = a s (m)/'K. The total physical 
cross section is obtained after incoherent addition of the contributions from 
the qq and gg production channels. 

Comparison of the partonic cross section at next-to-leading order with 
its lowest-order value reveals that the ratio becomes very large in the near- 
threshold region. Indeed, as r\ — ► 0, the "i-T-factor" at the partonic level K{r\) 
grows in proportion to a In 2 (rj) . The very large mass of the top quark, and the 
correspondingly small value of a notwithstanding, the large ratio K (77) makes 
it evident that the next-to-leading order result does not necessarily provide a 
reliable quantitative prediction of the top quark production cross section at the 
energy of the Tevatron collider. Analogous examples include the production of 
hadronic jets that carry large values of transverse momentum, the production 
of pairs of supersymmetric particles with large mass, and the pair-production 
of a fourth-generation quark, such as the postulated b' . 



2 Gluon Radiation and Resummation 

The origin of the large threshold enhancement may be traced to initial-state 
gluonic radiative corrections to the lowest-order channels. I remark that I am 
describing the calculation of the inclusive total cross section for the production 
of a top quark-antiquark pair, i.e., the total cross section for t + t + anything. 
The partonic subenergy threshold in question is the threshold for t + i + any 
number of gluons. This coincides with the threshold in the invariant mass of 
the t + i system for the lowest order subprocesses only. 

For i + j — ► t + i + g, the variable z is defined through the invariant 
(1 - z) = where k and p t are the four-vector momenta of the gluon and 

top quark. In the limit that z — > 1, the radiated gluon carries zero momentum. 
After cancellation of soft singularities and factorization of collinear singulari- 
ties in O(a^), there is a left-over integrable large logarithmic contribution to 
the partonic cross section associated with initial-state gluon radiation. This 
contribution is often expressed in terms of "plus" distributions. In O(al), it is 
proportional to a 3 In (1 — z). When integrated over the near-threshold region 
1 > z > 0, it provides an excellent approximation to the-full next-to-leading 
order physical cross section as a function of the top mass a. 

Although a fixed-order 0(a 4 ) calculation of ti pair production does not 
exist, universality of the form of initial-state soft gluon radiation may be in- 
voked, and the leading logarithmic structure at 0(a 4 ) may be appropriated 
from the next-to-next-to-leading order calculations of massive lepton-pair pro- 
duction (11), the Drell-Yan process. In the near-threshold region, the hard 
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kernel becomes 



2C 2 ln 4 (l-z)-^6 2 ln 3 (I-z) 



^° +1+2) (z,a) ~ l + 2aC y ln 2 (l-z) + a 2 

o 

(2) 

The coefScient 62 = (HCU — 2n/)/12; the number of flavors n/ = 5; Cqg = 
(Tp = 4/3; and C gg — Ca — 3. The leading logarithmic contributions in each 
order of perturbation theory are all positive in overall sign so that the leading 
logarithm threshold enhancement keeps building in magnitude at each fixed 
order of perturbation theory. 

The goal of gluon resummation is to sum the series in a n ln 2n (l — z) to all 
orders in a in order to obtain a more trustworthy prediction. This procedure 
has been studied extensively for the Drell-Yan process, and good agreement 
with data is achieved. In essentially all resummation procedures, the large log- 
arithmic contributions are exponentiated into a function of the QCD running 
coupling strength, itself evaluated at a variable momentum scale that is a mea- 
sure of the radiated gluon momentum. The set of purely leading monomials 
a™ln 2 "(l — z) exponentiates directly, with a evaluated at a fixed large scale 
fx = m, as may be appreciated from a glance at Eq. (|J). This simple result 
does not mean that a theory of resummation is redundant, even if only leading 
logarithms are to be resummed. Indeed, straightforward use of the exponen- 
tial of otlCij ln 2 (l — z) would lead to an exponentially divergent integral (and 
therefore cross section) since the coefficient of the logarithm is positive. The 
naive approach fails, and more sophisticated resummation approaches must be 
employed. 

Different methods of resummation differ in theoretically and phenomeno- 
logically important respects. Formally, if not explicitly in some approaches, 
an integral over the radiated gluon momentum z must be done over regions in 
which z — > 1. Therefore, one significant distinction among methods has to do 
with how the inevitable "non-perturbative" region is handled. 

The method of resummation employed in my work with Harry Contopana- 
is based on a perturbative truncation of principal-value (PV) resumma- 
This approach has an important technical advantage in that it does 
not depend on arbitrary infrared cutoffs. Because extra scales are absent, the 
method permits an evaluation of its regime of applicability, i.e., the region of 
the gluon radiation phase space where leading-logarithm resummation should 
be valid. We work in the MS factorization scheme. 

Factorization and evolution lead directly to exponentiation of the set of 
large threshold logarithms in moment (n) space in terms of an exponent E PV . 
The function E PV is finite, and lim n ^oo E pv (n, m 2 ) = —00. Therefore, the 
corresponding partonic cross section is finite as z — > 1 [n — > +00). 




3 



The function E pv includes both perturbative and non-perturbative con- 
tent. The non-perturbative content is not a prediction of perturbative QCD. 
Contopanagos and I choose to use the exponent only in the interval in moment 
space in which the perturbative content dominates. We derive a perturbative 
asymptotic representation of E(x 1 a(m)) that is valid in the moment-space 
interval 

1 < x = Inn < t = — — . (3) 
2ao 2 

The interval in Eq. (|3|) agrees with the intuitive definition of the perturbative 
region in which inverse-power contributions are unimportant: — !• 

The perturbative asymptotic representation is 

N(t) + 1 p+1 

E l0 (x, a) ~ Etj {x, a, N(t)) = 2C!y ^ of ■ (4) 

p=l j=0 

Here 

*J,P = -b^i-iY^Vc^-Ap ; (5) 

and r(l +z) = X^fcLo Cfe2 ^ wnere r i s the Euler gamma function. The number 
of perturbative terms N(t) in Eq. (^|) is obtained u by optimizing the asymp- 
totic approximation 



E(x,a)-E(x,a,N(t)) 



= minimum. Optimization works 



perfectly, with N(t) = 6 at m = 175 GeV. As long as n is in the interval of 
Eq. (||), all the members of the family in n are optimized at the same N(t), 
showing that the optimum number of perturbative terms is a function of t, i.e., 
of m only. 

Resummation is completed in a finite number of steps. When the running 
of the coupling strength a is included up to two loops, all monomials of the 
form a k ln fc+1 n, a k ln fc n are produced in the exponent of Eq. (0). We discard 
monomials a k ln fc n in the exponent because of the restricted leading-logarithm 
universality between tt production and massive lepton-pair production, the 
Drell-Yan process. 

The moment-space exponent that we use is the truncation 

N(t) + 1 

E ij (x,a,N)=2C ij a p s p xP +1 , (6) 

P =i 

with the coefficients s p = s p +i tP — b2~ 1 2 p / p(p+l). This expression contains no 
factorially-growing (renormalon) terms. One can also derive the perturbative 
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expressions, Eqs. ([?]), (^), and (|5|), without the principal- value prescription, 
although with less certitude 0. 

After inversion of the Mellin transform from moment space to the physi- 
cally relevant momentum space, the resummed hard kernel takes the form 

H%(z,a) = dxe Ei ^ x ' a ^y2Qj(x,a) . (7) 

The leading large threshold corrections are contained in the exponent Eij(x, a), 
a calculable polynomial in x. The functions {Qj(x,a)} arise from the ana- 
lytical inversion of the Mellin transform from moment space to momentum 
space. These functions are expressed in terms of successive derivatives of E. 
Each Qj contains j more powers of a than of x so that Eq. (Q) embodies 
a natural power-counting of threshold logarithms. However, only the lead- 
ing threshold corrections are universal. Final-state gluon radiation as well 
as initial-state/final-state interference effects produce subleading logarithmic 
contributions that differ for processes with different final states. Accordingly, 
there is no physical basis for accepting the validity of the particular subleading 
terms that appear in Eq. (Q). Among all {Qj} in Eq. (0), only the very leading 
one is universal, Qq, and it is the only one we retain. Hence, Eq. (^) can be 
integrated explicitly, and the resummed partonic cross sections become 

a* pert (r),m) = jf""" dze E ^ 1 ^M & '..(r 1 ,m,z). (8) 

The derivative <7^(r?,m,z) = d(a^\rj,m, z))/dz, and <7y is the lowest-order 
0(a 2 s ) partonic cross section expressed in terms of inelastic kinematic vari- 
ables. The lower limit of integration, z m i ni is fixed by kinematics. The upper 
limit, z m ax < 1, well specified within the context of our calculation, is estab- 
lished by the condition of consistency of leading-logarithm resummation. It 
is derived from the requirement that the value of all subleading contributions 
Qji j > 1 be negligible compared to the leading contribution Qq. The pres- 
ence of z ma x guarantees that the integration over the soft-gluon momentum 
is carried out only over a range in which poorly specified non-universal sub- 
leading terms would not contribute significantly even if retained. We cannot 
justify continuing the results of leading-logarithm resummation into the region 
1 ^ z > z max . 

To obtain the physical cross section, we insert the resummed expression 
Eq. (||) into Eq. (Q) and integrate over rj. Perturbative resummation probes the 
threshold down to rj > r/o = (1— z maa; )/2. Below this value, perturbation theory 
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pp->tt + X, Vs=1.8TeV 
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Figure 1: Inclusive total cross section for top quark production. The dashed curves show 
the upper and lower limits while the solid curve is our central prediction. CDF and DO data 
are shown. 

is not to be trusted. For m = 175 GeV, we determine that the perturbative 
regime is restricted to values of the subenergy greater than 1.22 GeV above 
the threshold (2m) in the qq channel and 8.64 GeV above threshold in the gg 
channel. The difference reflects the larger color factor in the gg case. The 
value 1.22 GeV is comparable to the decay width of the top quark, a natural 
definition of the perturbative boundary and by no means unphysically large. 

3 Physical cross section 

Other than the top mass, the only undetermined scales are the QCD factor- 
ization and renormalization scales. A common value /i is adopted for both. 
In Fig. 1, our total cross section for ii-production is shown as a function of 
top mass in pp collisions at y/S =1.8 TeV. The central value is obtained with 
the choice fi/m = 1, and the lower and upper limits are the maximum and 
minimum of the cross section in the range fi/m £ {0.5, 2}. At m — 175 GeV, 
the full width of this uncertainty band is about 10% . As is to be expected, less 
variation with /i is evident in the resummed cross section than in the next-to- 
leading order cross section. In estimating uncertainties, Contopanagos and I 
do not consider explicit variations of the non-perturbative boundary, expressed 



through z max . For a fixed to and /i, z max is obtained by enforcing dominance 
of the universal leading logarithmic terms over the subleading ones. Therefore, 
Zmax is derived and is not a source of uncertainty. At fixed to, the boundary 
necessarily varies as fi and thus a vary. 

Contopanagos and I calculate o u (jn — 175 GeV, y/S = 1.8 TcV) = 
5.52^Q 42 pb, in agreement with data EL This cross section is larger than 
the next-to-leading order value by about 9%. The top quark cross section 
increases quickly with the energy of the pp collider. We determine a (m = 
175 GeV, VS = 2 TeV) = 7.561°;*° pb. The central value rises to 22.4 pb at 
VS = 3 TeV and 46 pb at v 7 ^ = 4 TeV. _ 

Extending our calculation to larger values of m at \/S =1.8 TeV, we find 
that resummation in the principal qq channel produces enhancements over the 
next-to-leading order cross section of 21%, 26%, and 34%, respectively, for 
to = 500, 600, and 700 GeV. The reason for the increase of the enhancements 
with mass at fixed energy is that the threshold region becomes increasingly 
dominant. Since the qq channel also dominates in the production of hadronic 
jets at very large values of transverse momenta, we suggest that on the order 
of 20% of the excess cross section reported by the CDF collaboration!] may be 
accounted for by resummation. 



4 Other Methods of Resummation 

Two other groups have published calculations of the total cross section at m — 
175 GeV and y/S = 1.8 TeV: ct"(LSvN0) = 4.951°,;!°, Pb; and ct'^CMNTq) = 
4.75lQ g| pb. From a numerical point of view, ours and theirs all agree within 
their estimates of theoretical uncertainty. However, the resummation methods 
differ as do the methods for estimating uncertainties. Both the central value 
and the band of uncertainty of the LSvN predictions are sensitive to their 
arbitrary infrared cutoffs. To estimate theoretical uncertainty, Contopanagos 
and I use the standard /^-variation, whereas LSvN obtain theirs primarily from 
variations of their cutoffs. It is difficult to be certain of the central value and to 
evaluate theoretical uncertainties in a method that requires an undetermined 
infrared cutoff. 

The group of Catani, Mangano, Nason, and Trentadue (CMNT)a calculate 
a central value of the resummed cross section (also with /i/m — 1) that is 
less than 1% above the exact next-to-leading order value. Both they and 
we use the same universal leading-logarithm expression in moment space, but 
differences occur after the transformation to momentum space. The differences 
can be stated more explicitly if one examines the perturbative expansion of the 
resummed hard kernel HfAz : a). If, instead of restricting the resummation to 
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the universal leading logarithms, one uses the full content of H^(z, a), she or 
he would find an analytic expression that is equivalent to CMNT's numerical 
inversion, 



U% ~ 1 + 2aCij 



ln 2 (l-z) + 2 7jE ln(l-z) 



0(a 2 ). (9) 



In terms of this expansion, in our work we retain only the leading term ln 2 (l— z) 
at order a, but both this term and the non-universal subleading term 2je ln(l— 
z) are retained by CMNT. If this subleading term is discarded in Eq. ([)]), the 
residuals %/cr^ i0 defined by CMNT increase from 0.18% to 1.3% in the qq 
production channel and from 5.4% to 20.2% in the gg channel. After addition 
of the two channels, the total residual 5/a NLO grows from the negligible value 
of about 0.8% to the value 3.5%. While still smaller than the increase of 9% 
that we obtain, the increase of 3.5% vs. 0.8% shows the substantial influence 
of the subleading logarithmic terms in the CMNT results. 

Contopanagos and I judge that it is preferable to integrate over only the 
region of phase space in which the subleading term is suppressed numerically. 
Our reasons include the fact that the subleading term is not universal, is not the 
same as the subleading term in the exact 0(a 3 ) calculation, and can be changed 
if one elects to keep non-leading terms in moment space. The subleading term 
is negative and numerically very significant when it is integrated throughout 
phase space (i.e., into the region of z above our z max ). In our view, the results 
of a leading-logarithm resummation should not rely on subleading structures in 
any significant manner. The essence of our determination of the perturbative 
boundary z max is precisely that below z max subleading structures are also 
numerically subleading, whether or not these poorly substantiated subleading 
logarithms are included. _ 

In the remainder of this section I offer a more systematic analysis B of the 
role played in the CMNT approach by non-universal subleading logarithms 
and show in some detail how their method and results differ from ours. I treat 
expansions of the resummed momentum-space kernel up to two loops. The 
corresponding cross sections are integrable down to threshold, z max = 1 and 
rj = 0. However, the effects of the various classes of logarithms are pronounced 
if one continues the region of integration beyond our perturbative regime. 

In moment space, the exponent to two-loops is obtained from Eq. (^]): 

Ef^(x, a) = ga{s 2t ix 2 + s 1A x + s Qtl )+ga 2 (s 3t2 x 3 + s 2:2 x 2 + si t2 x + s Qt2 ), (10) 

with g — 2Cij and x = Inn. One can perform the analytical Mellin inversion 
directly, beginning with Eq. (|lO|). After a trivial integration, the results for 
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the one- and two-loop hard kernels are 



U {1) = x 2 z a{gs 2 ,i} + x z a{g{si A + 2cis 2: i)} , (11) 

and 

= x z a 2 {g 2 sl 1 /2} + x 3 z a 2 {gs 3 , 2 + g 2 (s 2 . lSl s + 2c lS ^)} 
+x 2 z a 2 {g(s 2 ,2 + 3cis 3! 2) + g 2 {s\ l /2 + 3ciSi,iS 2 ,i + s 2 ,is ,i + sj A [6c 2 - it 2 }} 
+x z a 2 {g(si i2 + 2cis 2 , 2 + s 3:2 [Gc 2 - ir 2 }) 
+5 2 («o,iSi,i + 2ciS ,lS2,l + ClS^! + s 2 4 s i,i[6c2 - 7T 2 ] 

+4 1 [12c 3 -27r 2 Cl ])} . (12) 

All the constants are defined in Eqs. (Q) and (||). Equation ( |il| ) includes a 
leading logarithmic term, cc^a, as well as a next-to-leading term, x z a. 

The question to be addressed is whether it is justified and meaningful to 
retain all of the terms in Eqs. ( |ll| ) and ( [l2| ) in the computation of the resummed 
cross section. The issue has to do with what one intends by resummation 
of leading logarithms. Contopanagos and I use the term leading logarithm 
resummation to denote the case in which the moment space exponent, Eq. ([To|), 
contains only the constants -Ell = {s p +i, p , 0}. This is also what is done in the 
CMNT method, and the exponent in moment space in their work is identical to 
that used for our predictions, Eq. (^). However, in contrast to our expression in 
momentum space, Eq. (^), the corresponding CMNT expression in momentum 
space includes the numerical equivalent of all terms in Eqs. ( p"l| ) and ( p^ ) that 
are proportional to s p+ i, p . 

If expressed analytically, CMNT's corresponding "LL" hard kernels are 



T~Lll = x l a 9 ~ x z a2g~f E , (13) 



and 



Uf h - x 4 z a 2 g 2 /2 + xla 2 {2gb 2 /3 - 2 lE g 2 } 
+x 2 z a 2 {-2gb 2lE + 5 2 [3 7 | - tt 2 /2]} 

+x z a 2 {2gb 2 [3 1 2 E - tt 2 /2]/3 + g 2 [ lE v 2 - 2 7 3 - 4C(3)]} , (14) 

where ((s) is the Riemann zeta function; C(3) = 1.2021)569. Evaluating the 
expressions numerically for the qq channel, one obtains El 

U [ l\ = x 2 z a x 2.66666 - x z a x 3.07848, (15) 
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and 



U { 1\ = x\a 2 x 3.55555 - x 3 z a 2 x 4.80189 

-x 2 z a 2 x 33.88456 - x z a 2 x 9.82479 . (16) 

Apart from the leading monomials that are the same as those in our approach, 
Eqs. ([[5]) and (|l^) include a series of subleading terms, each of which has a 
significant negative coefficient. In practice, these subleading terms suppress 
the effects of resummation essentially completely. One of the effects of this 
suppression is that the resummed partonic cross section is smaller than its 
next-to-leading order counterpart in the neighborhood of 77 = 0.1, a region 
in which the next-to-leading order partonic cross section takes on_its largest 
values. This point is illustrated in Fig. 3 of CMNT's second paperQ. 

Although the specific set of subleading terms in Eqs. (|l^) and ([h]) is gen- 
erated in the inversion of the Mellin transform, a case can be made that the 
terms are accidental. First, terms involving 7^ do not appear in the exact 
next-to-leading order calculation of the hard part, since they are removed in 
the specification of the MS factorization scheme. Therefore, the term pro- 
portional to je in Eq. (|l3| ) is suspect. Second, if the specific value of the 
subleading logarithm is extracted from the full 0(a 3 ) next-to-leading order 
calculation, one finds u x z a(2g — 41/6) instead of the term ~x z a2gjE- In- 
stead of the numerical coefficient 3.07848 in Eq. (|l5|), one finds the smaller 
value 1.5 if the subleading logarithm of the exact 0(a 3 ) calculation is used. 
Thus, not only is the O(a) subleading term retained in the CMNT approach 
different from that of the exact calculation, it is numerically about twice as 
large. Third, the results of a LL resummation should not rely on the sub- 
leading structures in any significant way. However, in the CMNT approach, 



Eq. (13), which is the one-loop projection of their resummed prediction, repro- 
duces only 1/3 of the exact 0(a 3 ) enhancement, the other 2/3 being cancelled 
by the second (non- universal) term of Eq. ( |l3| ) . Although the goal is to resum 
the threshold corrections responsible for the large enhancement of the cross 
section at next-to-leading order, the CMNT method does not reproduce most 
of this enhancement. 

Addressing questions associated with the^s terms, CMNT examine a type 
of NLL resummation in their second paper □. In this NNL resummation, the 
{■ s p+i,pi s p,p\ terms are retained in the exponent, Eq. (|To|). The corresponding 
hard kernels become 

n< NLL = x t a 9, (17) 
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and 



ri§LL = xta 2 gy2W z a 2 2gb 2 /3-xla 2 g 2 [ 1 l+7r 2 /2}-x z a 2 {gb 2 [2 1 2 E +7r 2 /3}+gHa3)} 

(18) 

Equation (|17|) is identical to the one-loop projection of our hard kernel. On 
the other hand, our two-loop projection contains only the first two terms of 
Eq. (|l8|). The term proportional to x^a 2 is present in our case, along with 
the leading term proportional to x z a 2 , because it comes from the leading 
logarithms in the exponent E(n), through two-loop running of the coupling 
strength. In contrast to Eq. (H), Eq. ( jig ) relegates the influence of the am- 
biguous constant coefficients to lower powers of x z (but with larger negative 
coefficients). In the amended scheme, the unphysical 75 terms are still present 
in the two-loop result, Eq. (|T^), along with ir 2 and £(3) terms that may be 
expected but whose coefficients have no welLdefined physical origin. Recast in 
numerical form, Eqs. (O) and dlq) becomeO 



Hnll = x t a x 2.66666, (19) 



and 



T-L% ] LL = x\a 2 x 3.55555 + 2^a 2 x 3.40739- z 2 a 2 x 37.46119-x z a 2 x 54.41253 . 

(20) 

There is a significant difference between the coefficients of all but the very 
leading power of x z in Eqs. Jl5| ) and ( [l6| ) with respect to those in Eqs. (19) 



and d20| ) , and the numerical coefficients grow in magnitude as the power of x 
decreases. 

Using their NLL amendment, CMNT find that the central value of their 
resummed cross section exceeds the next-to-leading order result by 3.5% (both 
qq and gg channels added) . This increase is about 4 times larger than the cen- 
tral value of the increase obtained in their first method, closer to our increase of 
about 9%. The reason for the significant change of the increase resides with the 



subleading structures, viz., in the differences between the LL version Eqs. (15) 
and ( |l6|) and the NLL version Eqs. ( |l9| ) and (^0|). The subleading terms at 
two- loops cause a total suppression of the two- loop contribution (in fact, that 
contribution is negative), if one integrates all the way into what we call the 
non-perturbative regime. This suppression explains why an enhancement of 
only 3.5% is obtained in the amended method, rather than our 9%. 

CMNT argue that retention of their subleading terms in momentum space 
is important for "energy conservation" . By this statement, they mean that one 
begins the formulation of resummation with an expression in momentum space 
containing a delta function representing conservation of the fractional partonic 
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momenta. In moment space, this delta function subsequently unconvolves the 
resummation. Therefore, when one inverts the Mellin transform to return to 
momentum space, the full set of logarithms generated by this inversion are 
required by the original energy conservation. This line of reasoning would be 
compelling if the complete exponent E(n) in moment space were known exactly, 
i.e., if the resummation in moment space were exact in representing the cross 
section to all orders. However, the exponent is truncated in all approaches, 
and knowledge of the logarithms it resums reliably is limited both in moment 
and in momentum space. Hence, the set of logarithms produced by the Mellin 
inversion in momentum space should also be restricted. In our approach en- 
ergy conservation is obeyed in momentum space consistently with the class of 
logarithms resummed. On the other hand, in the CMNT method, knowledge is 
presumed of all logarithms generated from the Mellin inversion, despite the fact 
that the truncation in moment space makes energy conservation a constraint 
restricted to the class of logarithms that is resummable, i.e., a constraint re- 
stricted by the truncation of the exponent E(n). The two approaches would 
be equivalent provided a constraint be in place on the effects of subleading 
logarithms. This constraint is precisely our restriction z max < 1, but no such 
constraint is furnished by CMNT. For this reason their results change signif- 
icantly if one set of the logarithms generated in momentum space is adopted 
as "the set corresponding to energy conservation" , and then compared with 
another set, produced by a different truncation of E(n). 

The essence of our determination of the perturbative regime, z max < 1, 
is precisely that, in this regime, subleading structures are also numerically 
subleading, whether or not the classes of subleading logarithms coming from 
different truncation of the master formula for the resummed hard kernel are 
included. The results presented in Fig. 11 of our second paperl3, show that if 
we alter our resummed hard kernel to account for subleading structures but 
still stay within our perturbative regime, the resulting cross section is reduced 
by about 4%, within our band of perturbative uncertainty. 

A criticism □ is that of putative "spurious factorial growth" of our re- 
summed cross section, above and beyond the infrared renormalons that are 
eliminated from our approach. The issue, as demonstrated in Eq. (29) of our 
second paper a, can be addressed most easily if one substitutes any monomial 
appearing in Eq. (|T^), symbolically a m c(l, m) In' x z , into Eq. (||) and integrates 
over z: 

/i l 
dz\n l x z =a m c{l,m)(l-z min )l\J2 l v j {l/{l-Zrnin)) • (21) 

For the purposes of this demonstration we set <r- ■ = 1. The coefficients c(l, m) 
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can be read directly from Eq. (|l2|). For the leading logarithmic terms, 

c(2m,m) oc 1/ml, (22) 

where this factorial comes directly from exponentiation. After the integration 
over the entire z-range, the power of the logarithm in x z becomes a factorial 
multiplicative factor, l\. The presence of l\ follows directly from the existence 
of the powers of In x z that are present explicitly in the finite-order result in 
pQCD and is therefore inevitable. If this exercise is repeated, but with the 
range of integration in Eq. (|2l]) constrained to our perturbative regime, one 



obtains the difference between the right-hand-side of Eq. (21) and a similar 
expression containing z max . The result is numerically smaller, but both of the 
pieces are multiplied by ^!. 

The factorial coefficient V. is not the most important source of enhance- 
ment. For the leading logarithms at two-loop order, I = 2m — 4, and the 
overall combinatorial coefficient from Eqs. ( pl| ) and ( p2| ) is (2m)!/m! = 12. 
For comparison, at representative values of r\ near threshold, r\ = 0.1 and 
0.01, the sum of logarithmic terms in Eq. ( pl| ) provides factors of 16.1 and 
314.3, respectively. Similarly, the (multiplicative) color factors at this order of 
perturbation theory are (2CV,) 2 = 7.1 and 36 for the qq and gg channels, re- 
spectively. All of these features are connected to the way threshold logarithmic 
contributions appear in finite-order pQCD and how they signal the presence 
of the noiuperturbative regime. Thus, preoccupation with the l\ factor seems 
misplacedla. 

"Absence of factorial growth" is based on the use by CMNT of Eq. ( |l6| ) for 
their main predictions, an expression that contains non-universal subleading 
logarithms, all with significant negative coefficients. Mathematically, factorial 



growth is present for each of the powers of the logarithm in Eq. (21), since these 
monomials are linearly independent. Absence of factorial growth based on a 
numerical cancellation between various classes of logarithms, most of them with 
physically unsubstantiated coefficients, appears to us to be an incorrect use of 
terminology. In the CMNT approach the effects of resummation are suppressed 
by a series of subleading logarithms with large negative coefficients. If there is 
no physical basis for preference of Eqs. (|l3|) and ( [hi]) over Eqs. (|lj) and (^8|), 
as CMNT appear to suggest, then the difference in the resulting cross sections 
can be interpreted as a measure of theoretical uncertainty. This interpretation 
would not justify firm conclusions of a minimal 0.8% increment in the physical 
cross section based on the choice of Eqs. ([b^) and (|l4|). 

The CMNT value for the inclusive top quark cross section at m = 175 GeV 
and ^/S = 1.8 TeV, including theoretical uncertainty, lies within our uncer- 
tainty band. Therefore, the numerical differences between our results for top 
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quark production at the Tevatron have little practical significance. However, 
there are important differences of principle in our treatment of subleading con- 
tributions that will have more significant consequences for predictions in other 
processes or at other values of top mass and/or at other energies, particularly 
in reactions dominated by gg subprocesses. 

5 Discussion and Conclusions 

The advantages of the perturbative resummation method i are that there are 
no arbitrary infrared cutoffs and there is a well-defined region of applicability 
where subleading logarithmic terms are suppressed. When evaluated for top 
quark production at yS — 1.8 TeV, our resummed cross sections are about 
9% above the next-to-leading order cross sections computed with the same 
parton distributions. The renormalization/factorization scale dependence of 
our cross section is fairly flat, resulting in a 9 — 10% theoretical uncertainty. 
Our perturbative boundary of 1.22 GeV above the threshold in the dominant 
qq channel is comparable to the hadronic width of the top quark, a natural 
definition of the perturbative boundary. 

Our estimated theoretical uncertainty of 9 — 10% is associated with \x 
variation. An entirely different procedure to estimate the overall theoretical 
uncertainty is to compare our enhancement n of the cross section above the 
next-to-leading order value to that of CMNTq, again yielding about 10%. An 
interesting question is whether theory can aspire to an accuracy of better than 
10% for the calculation of the top quark cross section. To this end, a mastery of 
subleading logarithms would be desirable, perhaps requiring a formidable com- 
plete calculation at next-to-next-to-leading order of heavy quark production, 
to establish the possible pattern of subleading logarithms, and resummation 
of both leading and subleading logarithms. An analysis in moment space of 
the issues involved in resummation of next-to-leading logarithms for heavy 
quark production is presented by Kidonakis and StcrmanH. Inversion of the 
resummed moments to the physically relevant momentum space requires con- 
siderable work. Full implementation of the resummation of next-to-lcading log- 
arithms would reduce the difference somewhat between our results and those of 
CMNT and move the debate to the level of next-to-next-to-leading logarithms. 

Our prediction falls within the relatively large experimental uncertainties. 
If a cross section significantly different from ours is measured in future ex- 
periments at the Tevatron with greater statistical precision, we would look 
for explanations in effects beyond QCD perturbation theory. These explana- 
tions might include unexpectedly substantial non-perturbative effects or new 
production mechanisms. An examination of the distribution in r\ might be 
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revealing. 

The all-orders summation of large logarithmic terms, that are important 
in the near-threshold region of small values of the scaled partonic subenergy, 
r\ — > 0, was described here for the specific case of top quark production at the 
Fermilab Tevatron collider. Other processes for which threshold resummation 
will also be pertinent include the production of hadronic jets that carry large 
values of transverse momentum and the production of pairs of supersymmetric 
particles with large mass. 
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